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1. Introduction 

The very accurate experimental measurements of the mass differences between the heavy and 
light B° s and B d mass eigenstates, AM S [jl]] and AM d [§], that describe the - fi° and B° d - B d mix- 
ings respectively, make improving the theoretical study of these quantities crucial. In the standard 
model (SM), the mass difference is given by [^] 

AM s{d) \ theor . = t^K w v tb \ 2 n2s (x,)M B ^j$^6 B<i) , (l.i) 

where x t = mj/M^,, r}^ * s a perturbative QCD correction factor, So(x t ) is the Inami-Lim function 
and the products fB s{d J*B s{f[] parametrize the hadronic matrix elements in the effective theory with 
fB s{d) the decay constants and Bg s(d) the (renormalization group invariant) bag parameters. The 
hadronic matrix elements can be calculated in lattice QCD. Our current knowledge of them limits 
the accuracy with which the CKM matrix elements appearing in Eqn. ( |1.1[ ) can be determined from 
the experimental measurements of AM s ( d y The goal of our project is to calculate all the hadronic 
matrix elements which are relevant for the mass and width differences in the B®, d ^ systems in 
unquenched lattice QCD at the few percent level. 

Many of the uncertainties that affect the theoretical calculation of the decay constants and 
bag parameters cancel totally or partially if one takes the ratio E, 2 = Ib^bJ fs d ^B d - Hence, this 
ratio and therefore the combination of CKM matrix elements related to it by Eqn. dl.l| ) can be 
determined with a significantly smaller error than the individual matrix elements. This is a crucial 
ingredient in the unitarity triangle analysis. In these proceedings we report our preliminary results 
for the determination of i§ , as well as for the quantities /J B B . 

Other work on this subject using 2+1 lattice QCD methods can be found in [Q]. 

2. Operators, actions and matching calculation 

The whole set of operators whose matrix elements are needed to determine the B®, ^ mixing 
parameters are 

Q s[d) = [b%(i-75)s i (d i )} [£V(i_ 75 y (</;')] , 
ef ) = [b\\ - y 5 y(^')] [^'(i - K)sW] , 

Q f ) = [tf (1 - y 5 >V)] [#(1 " 75)s i (d i )] , (2-1) 

where i,j are color indices. In these proceedings we focus on the results for the first two pairs of 
operators, enough to determine AM S ^, and leave the study of the third pair, needed for an improved 
determination of Ar^), for a forthcoming publication 

We use the Fermilab action ^ for the b valence quarks and the Asqtad action for the light 
sea and valence quarks, u, d and s. The Fermilab action has errors starting at 0(a s Agco/M) and 
0((Aq CD /M) 2 ), while the errors of the Asqtad action are 0(a s a 2 ,a 4 ). 

The products /J Bg in Eqn. ( |1.1[ ) parametrize the matrix elements by 

(B° s \^\B° s ) m (n) = \Ml^^{ji) . (2.2) 
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The lattice matrix elements {^^\Q \^^) lat determine /f (rf) ^B l(rf) at tree level. Beyond tree- 
level, the operators Q s ( d \ mix with both on the lattice and in the continuum. Including 
one-loop corrections, the renormalized matrix element is given by 



^ _ (g*0)A«( M ) = [\ + a s -p L ^,m b )}(Q^) ht (a) + a s -p LS ^,m h )(Qf >) lat (a). (2.3) 



The O ( -jjjpH improvement is implemented by a rotation of the b quark as explained in [g], so 
the perturbative matching errors start at 0(aj,a s A/aM). The matching coefficients pn and pis 
are the differences between the continuum MS and lattice renormalization coefficients calculated at 
one-loop order. We have calculated these coefficients for the same choice of lattice actions as used 
in the numerical simulations. We have checked that our results have the correct infrared behavior, 
that they are correct in the massless limit, and that they are gauge invariant. However, our results for 
the matching coefficients are still preliminary, because not all diagrams have been independently 
checked. 



The optimal value of the strong coupling constant to be used in Eqn. (2.3) is the V-scheme 
coupling (Xy{q*) [||] where q* is defined in In the absence of a calculation of q* for the 

specific processes we are studying, we choose q* =2/ a, very close to the q* calculated for heavy- 
light currents. The specific values for a s we use are given in Table [j] 

3. Simulation details 

The matrix elements needed to determine both fl Bn and Be , are extracted from the follow- 
ing three-point and two-point functions 

Co(^i^2)=£<fe(^?i)r5?(^fi)|O(0)|fc(y,/ 2 )75 9 (y,f 2 )) ; 
C z (t)=£(b(x,t)Y 5 q(x,t)q(0)Y5b(0)), C A4 (t) = J £(Hx,t)Y 7 5 q(x,t)q(0)75b(0)), (3.1) 



where the operator O is any Q s ^ or Qg defined in Eqn. (2A ). The B meson operators are smeared 
at the sink with a 1 S onium wavefunction. All the correlation functions in Eqn. ( 3. 1 ) are calculated 
using the open meson propagator method described in []TC|]. 

We have performed these calculations on the MILC coarse lattices (a = 0.12 fm) with 2+1 
sea quarks and for three different sea light quark masses. The strange sea quark mass is always 
set to 0.050. The light sea masses, mi = mf^ ht , number of configurations and other simulation 
details are collected in Table ffl. The mass of the bottom quark is fixed to its physical value, while 
for each sea quark mass we determine the different matrix elements for six different values of 
the light valence quark mass in a generic meson B° q , m q = 0.0415,0.03,0.020,0.010,0.007,0.005. 
We use m q = 0.0415 for the valence strange mass in our simulations. It is close to the physical 



strange quark mass, m p s hys ' = 0.036 JTT|]. The matrix elements of the operators are extracted from 



simultaneous fits of three-point and two-point functions using Bayesian statistics. 
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mi/mf a 


Volume 


confs 


fl-^GeV) 


a s = a v (2/a) 


ZY sources 


0.020/0.050 


20 3 x 64 


460 


1.605(29) 


0.31 


4 


0.010/0.050 


20 3 x 64 


590 


1.596(30) 


0.31 


4 


0.007/0.050 


20 3 x 64 


890 


1.622(32) 


0.31 


4 



Table 1: Simulation parameters and a s used in the matching with the continuum, w/ is the light sea quark 
mass. 



4. Results 



The results in Fig. [T] show fg q J B^ s (nib) in lattice units as a function of the light valence mass 
am q . The errors shown are statistical errors only; the analysis of the systematic errors is not yet 
complete. 
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Figure 1: fe q y (»%) in lattice units. The different symbols and colors correspond to different values of 
the sea light quark masses, m;. 

The statistical errors range between 1 — 3%. Some conclusions can be already extracted from 
this plot. The light sea quark mass dependence of fg q \J^B q is small compared to the statistical 
errors. The dependence on the light valence quark mass, however, is noticeable within statistics. In 

order to get a value for fg j J 'b^ s {nib), since the s valence quark mass we are using is slightly larger 
than the physical one, we need an interpolation in the s valence quark mass together with a chiral 
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extrapolation to the physical sea quark masses. To determine fy d y B^(mb) we need extrapolations 
in both the d valence quark mass and the sea quark masses. This is in progress. 

In Fig. ^| we plot the ratio £ = fg s \JB#J fB d \/^B d as a function of the d valence quark mass 
in the denominator. Again, our results are preliminary for the same reasons as mentioned before 
and the errors are only statistical. Most of the systematic errors cancel in the ratio, but not those 
associated with the chiral extrapolation in the light valence quark mass. 
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Figure 2: | as a function of the valence d mass for three different values of the light sea quark masses. 

4.1 Chiral extrapolation 

The continuum chiral expansion of the hadronic matrix element (B q \Q\B q ^ at NLO in (par- 
tially quenched) heavy meson chiral perturbation theory (HMChPT) is given by Jl2| ] 

(B q \Q\B q ) = p(l+w(T q + W q + S q ))+c m q + c 1 {m u + m D +m s ) (4.1) 

where m\j ,mj),ms are the sea quark masses and m q the light valence quark mass. j8, w, Co and 
c\ are low energy constants (LECs) to be determined from the fits. The functions T q , W q , and S q 
contain the chiral logs and correspond to tadpole-, wave-function, and sunset-type contributions, 
respectively. 

The effects of 0(a 2 ) taste changing interactions can be included in Eqn. ( fi~I| ) using staggered 
chiral perturbation theory (SChPT). In that case, the chiral log functions are modified to depend 
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statistics 
scale(a~ 1 ) 
Higher order matching 
Heavy quark discret. 
Light quark discret. + ^PT fits 



fB q y/W q % 

1- 3 1-2 
0.9 

~ 4.5 cancel to a large extent 

2- 3 <0.5 

Work in progress 



Table 2: Error budget for /s \/Bb„ and % in percent. 



on the masses of the different taste multiplets. Explicit expressions from SChPT for heavy-light 



bilinear quantities can be found in [0, 14]. Similar terms are expected to contribute to our four- 



quark operators. The modified chiral logs contain other fixed LEC's, most of which are already 



determined to a high degree of certainty [14]. The logs also contain constants from heavy quark 
effective theory, in particular the mass splitting between the vector and pseudoscalar heavy mesons, 
A* = Mb * —Mb, and the mass splittings between the pseudoscalar heavy mesons containing dif- 
ferent valence and sea light quarks, 8 qr = Mg r —Mb . These HQET constants can be determined 
directly from the two-point function fits and used as input into the chiral fits, with the experimental 
values then used in the extrapolation. 



We are still in the process of determining the exact SChPT form of Eqn. (|4. 1|). Once the 



functional SChPT form of Eqn. ( |4.1[ ) is completely determined, we plan to use it to simultaneously 
fit it to our lattice data points for all sea and valence quark masses, and to determine the unknown 
LEC's in the process. For the systematic error analysis, we plan to study the effects of changing the 
SChPT form, for example by adding NNLO analytic terms, and the effects of allowing the more 
poorly known fixed parameters to vary. Our physical results will then be obtained by turning the 
taste- violations off, and extrapolating (interpolating) to the physical light (strange) sea and valence 
quark masses. 

We will quote results for t, , fs s \JBb s , fs d y/BB d , Bb s , and Bs d when this step is completed. We 
expect light quark discretization effects to be an important source of uncertainty until we calculate 
the three-point correlators at several lattice spacings and use these in the chiral fits. 

5. Summary and future work 

We have presented preliminary results for fs q ^Bs q for six different values of m q as well as 
for the ratio E, with five different values of m v d a,ence . Our analysis on three ensembles with different 
sea light quark masses gives statistical errors between 1 — 3% for fs q \jB~B~ q and 1 — 2% for the ratio 
t, . The systematic error analysis is in progress, see Table |2[ 

Two important sources of error in the calculation of fs q \JBs q , the matching uncertainties and 
heavy quark discretization errors, are expected to cancel to a large extent when taking the ratio. 
We have already checked that the difference between the tree level and the one-loop results for ^ 
with our preliminary results for the renormalization constants is less than 0.5%. The higher order 
matching errors in Table |2] have been naively estimated as being 0(l x af) for the coarse lattice. 
Heavy quark discretization effects in the table are estimated by power counting [|15|] . 
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We are in the process of generating lattice data on the coarse lattice with a smaller light sea 
quark mass, mi = 0.005, which will further constrain the chiral extrapolations. Better fitting ap- 
proaches and different smearings that could reduce statistical errors further are also being inves- 
tigated. Results from these improvements will be presented in a future publication We also 
plan to present results for the decay width differences Ar 5 and Ar^, for which we have already 
calculated all the hadronic matrix elements needed. Finally, we plan to improve this analysis by 
repeating this calculation at other lattice spacings to study the discretization errors in detail. Simu- 
lations on finer lattices will reduce both discretization and perturbative matching errors. 
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